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ABSTRACT 

Unsteady aerodynamic forces are calculated by the XTRAN2L finite 
difference program which solves the complete two-dimensional unsteady tran- 
sonic small perturbation equation. The unsteady forces are obtained using a 
pulse-transfer function technique which assumes the flow field behaves in a 
locally-1 inear fashion about a mean condition. Forces are calculated for a 
linear flat plate using the default grids from the LTRAN2-NLR, LTRAN2-HI, and 
XTRAN3S programs. The forces are compared to the exact theoretical values for 
flat plate, and grid-generated boundary and internal numerical reflections are 
observed to cause significant errors in the unsteady airloads. Grids are pre- 
sented that alleviate the reflections while reducing computational time up to 
fifty-three percent and program size up to twenty-eight percent. Forces are 
presented for a six percent thick parabolic arc airfoil which demonstrate that 
the transform technique may be successfully applied to non-linear transonic 
flows. 


NOMENCLATURE 

airfoil chord 

unsteady lift force coefficient, per radian of pitch angle 
unsteady pitching moment coefficient, per radian of pitch angle 

P-P«, 

pressure coefficient, 

^00 

airfoil contour 
reduced frequency, 

* 00 

freestream Mach number 
local static pressure 
freestream static pressure 
freestream dynamic pressure 
time, seconds 
freestream velocity 

orthogonal coordinates in physical space 

X z 

non-dimensional orthogonal coordinates in physical space (— — ) 

angle of attack radians 
ratio of specific heats 
2 - {2 -y)M» 2 
airfoil thickness ratio 

orthogonal coordinates in computational space 
scaled time, 


c 

cm„ 


f 

k 

M« 

P 

poo 

Qoo 

t 

Vo, 

X,z 

x,z 

a 

Y 

Y* 

5 

5.n 


T 



At scaled time step 

((> perturbation velocity potential 

w frequency of oscillation, radian/second 

Subscript 

0 initial or steady-state value 

[ ] denotes jump in quantity across a discontinuity 


INTRODUCTION 

Considerable effort is underway to develop finite difference computer 
codes for calculating transonic unsteady aerodynamics for flutter and other 
aeroelastic analyses. Difference equations are solved for a finite number of 
discrete points in the flow field using the techniaues that have enjoyed 
considerable success in steady transonic flow analyses. The distribution of 
grid points for steady flows is a topic of active research. ^ Current 
practice is to map the physical domain to a finite computational region using 
smoothly varying stretching functions. In the calculation of unsteady flows 
an added concern is the dynamic behavior of the computed solution since 
spurious reflections due to rapid changes in grid spacing and from the outer 
edges of the computational domain are possible. In steady flows these effects 
are suppressed by iterating solutions until all dynamic phenomenon have 
subsided. Thus special emphasis must be given to the development and 
evaluation of finite difference grids for unsteady problems. 

Several investigations have calculated unsteady transonic flows with 
finite difference techniques. ^"2^ Some have used the higher level flow 
equations such as the Navier-Stokes equations^, the Euler equations^"^ and 
the full potential equation. 6-11 In these codes the grid is generally 
curvilinear, moving with the airfoil motion and ouite complex in nature. In 
these works some variations of the internal parameters are made to verify the 
results but little documentation of the influence of the grid variations is 
given. Considerable work has been done using the low freouency versions of 
the transonic small oerturbation (TSP) eauationl^"!^ including extensive 
flutter calculations. 16“17 Recent efforts have been directed at solving the 
complete TSP equation^ which is applicable to higher reduced 
frequencies. The grid system for the TSP equations is considerably simplified 
as the boundary conditions are applied on a mean line representing the airfoil 
and wake, and the grid remains stationary. Normally a rectilinear grid is 
used with various spacing schemes to enhance the treatment of local phenomena. 

There are some documented efforts that have shown significant effects of 
the grid on unsteady flow problems. 21 In ref. 11, it was demonstrated 
that if the grid was stretched too rapidly, such that the flow field was 
inadequately resolved, spurious effects can occur in the time domain 
solution. The investigation of ref. 20 indicated that such effects could be 
alleviated for the TSP equation with a more gentle stretching of the grid. 
During the development of the freouency domain perturbation method of refs. 
21-22, significant difficulty was encountered until it was noted that an 
inadequate number of field points were being used to describe the wavelengths 
involved. Increasing the point density corrected the problem. 
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In the initial application of the three-dimensional TSP program 
XTRAN3 s 23 at NASA Langley Research Center, indications of grid generated 
anomalies were observed. To permit economical assessment of the grid 
questions, a corresponding two-dimensional program, XTRAN2L, was used. The 
program XTRAN2L is a version of the HYTRAN2 codel® modified to solve the 
complete TSP equation. The purpose of this paper is to describe both the 
results and the methodology used in assessing the computational grids. The 
approach taken is basically one of performing numerical experiments. The 
linearized TSP equation applied to a flat plate airfoil is used for the test 
case in order to permit comparisons with exact results. The airfoil is given 
a small prescribed pulse in angle of attack and the aerodynamic transients are 
calculated. A transfer function analysis using fast Fourier transforms is 
then used to obtain a detailed description of the aerodynamic forces in the 
frequency domain. This paper presents the results for various computational 
grids and also illustrates their use for a parabolic arc airfoil in transonic 
flow. 


ANALYSIS TECHNIQUE AND APPLICATION 

Unsteady Transonic Small Perturbation Equation 

All calculations were obtained using the code XTRAN2L which solves the 
complete unsteady transonic small perturbation (TSP) potential equation 


? ? 2 
k M 2 kM„ 


-M 

7 ^ 




( 1 ) 


where t* = 2 - (2-y)Mo92,. which is based on the scaling used in the 
LTRAN2-NLR program. The airfoil flow tangency and trailing edge 

conditions are 'applied on the n = 0 line and, in the small perturbation 
approximation, become 


= f^ + kf^ ;n=+ 0 , 0£5£l 

( 2 a) 

[♦ 5 ] + =0 ; n = 0 , 5 = 1 

( 2 b) 


The wake, which is represented as a slit downstream of the airfoil trailing 
edge, has the boundary condition 


[♦ 5 ] + kUx] = 0 ; n = 0, ? > 1 (3) 

Numerical solutions of Eq. (1) were obtained using the alternating- 
direction implicit algorithm of Rizzetta and Chin. 18 The Rizzetta-Chin 

algorithm is similar to that used in the LTRAN2 codeia with the addition 
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of a three-time-level representation of the 


term on the n-sweep. 

The code XTRAN2L was developed by incorporating, the Rizzetta-Chin algorithm 
into the low /moderate frequency code HYTRAN2.16 The HYTRAN2 program is a 
modified version of the LTRAN2-NLR code which was derived from the LTRAN2 
program. The TSP equation and boundary conditions in the LTRAN2 program are 
low frequency approximations and do not include the time derivative terms. 
The HYTRAN2 code (and LTRAN2-NLR) is a moderate frequency version and includes 
the time derivative terms in the boundary conditions. Since XTRAN2L solves 
the complete TSP equation, it has no restrictions on the allowable values of k 
(LTRAN2 begins to fail at k « 0.075, and HYTRAN2 is valid for k < 0.4). 


Engquist and Majda^S developed far-field radiation boundary conditions 


for the low frequency equation (without the 




term) which were 


incorporated into LTRAN2 by Kwak.26 Those boundary conditions reduced the 
disturbances reflecting from the computational boundaries. Since the boundary 
conditions of ref. 25 are applicable only to the low freouency equation, 
far-field radiation conditions consistent with Eq. (1) have been developed and 
implemented into XTRAN2L. These conditions allow the use of smaller grids 
while significantly reducing reflections from the grid boundaries. Far-field 
radiation boundary conditions were used in analyzing the XTRAN2L default grid 
while reflecting far-field boundary conditions described in ref. 14 were used 
in analyzing all the other grids. 


Pulse Technique 


In order to determine the accuracy of results obtained using a particular 
grid, the linear unsteady aerodynamic forces on a flat plate are calculated 
for a wide range of reduced frequencies. The forces are compared to the known 
exact theoretical values to determine the grid characteristics. Typically, 
unsteady aerodynamic forces are determined by calculating several cycles of 
forced harmonic oscillations with the last cycle providing the estimate of the 
forces. Alternatively, the estimates may be obtained from Fourier transforma- 
tion of the response to a step change in a given mode of motion. The 
harmonic oscillation technioue is costly since calculations may be required 
for 5-10 reduced frequencies. In contrast, the indicial response approach can 
be inaccurate for finite difference solutions of Eq. (1). The discontinuous 
initial condition caused by a step function can only be roughly approximated 
and thus an error in the calculation of the airfoil loads is introduced and 
the resulting forces can be inaccurate. It was suggested in ref. 27 that the 
use of a smooth pulse would alleviate this problem. A different pulse from 
that of ref. 27 is used in this study. It is a pulse in pitch about the 
quarter-chord defined by 

n e-(^“‘fc)^ 

“ * 360 

where 

Tc = 57.5 at 

This gives a pulse with a maximum amplitude of one-half degree with a 
smooth transition from a steady-state condition to the pulse. The calculation 
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is carried out for a total of 1024 time steps to assure a return of the flow 
field to the steady-state condition. 

Fast Fourier transforms (FFTs) of the lift coefficient, the moment 
coefficient about the quarter-chord, and the angle of attack time histories 
are then calculated. The lift and moment coefficient FFTs are divided by the 
o FFT to obtain the frequency response functions for and Cm^^ as 

shown in Fig. 1. For most cases the unsteady lift and moment show the same 
grid characteristics, and only the calculated lift is shown. The pulse was 
designed to have sufficient magnitude to obtain reasonable results for 
frequencies up to k = 2.0. The use of the pulse and transfer function 
technique gives considerable detail in the frequency domain with a significant 
reduction in cost over calculating discrete oscillations. 

The linear analysis is performed for a flat plate at M« = 0.850, 
At = it/12.80, ao = 0. rad. and V«,/c = 100.0. The resulting forces are 
compared to the exact forces calculated by a two-dimensional subsonic kernel 
function aerodynamics program described in ref. 28. The comparison gives an 
excellent indication of the effect of the grid for a wide range of reduced 
frequencies. The analysis is also performed on a 6 % thick parabolic arc 
airfoil at the same conditions as the flat plate to demonstrate the non-linear 
transonic effects. For lifting conditions, the steady-state and unsteady 
algorithms converge to slightly different solutions, so before the pulse is 
applied the airfoil is held fixed and the unsteady solution is calculated for 
1024 time steps to obtain a consistent initial condition. This ensures that 
the lift and moment after the pulse return to the initial values and 
eliminates low frequency errors which might be introduced in the FFT results. 


Grids Analyzed 


Five computational grids are evaluated using the pulse-transfer function 
technioue. Three are the default grids from currently used TSP finite 
difference programs; the LTRAN2-NLR prooram,^^ the LTRAN2-HI program, 
and the x,z grid of the XTRAN3S program. 23 In addition, two grids designed 
to improve the accuracy of the calculated forces are analyzed. 

The default grid of the LTRAN2-NLR program consists of 99 x 79 points in 
x,z, and the grid coordinates for the flat plate case are listed in Table 1. 
Figure 2 illustrates the LTRAN2-NLR grid in the region around the flat plate 
with a typical airfoil drawn to illustrate the airfoil location. The physical 
grid for the flat plate case extends from -1034c to 856c in x and +3860c in z 
with 33 grid points on the airfoil. Since the far-field boundary conditions 
of the LfRAN2-NLR program are reflecting, the physical mesh extent is several 
thousand chordlengths in both directions. The TSP equation is solved in a 
transformed coordinate system, so the physical grid is mapped into 
computational coordinates. In LTRAN2-NLR the computational space C,n grid is 
identical for all cases and related to the physical x,z grid by 


The physical x grid is fixed and the z grid varies with Mach number, thickness 
ratio, and ratio of specific heats. For the flat plate airfoil a value of 
6 = 0.005 is used. 
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TABLE 1. LTRAN2-NLR DEFAULT GRID FOR M« = 0.850, 6 = 0.005 


INDEX 

X 

2 

1 

-1033.53047 

-3860.13477 

2 

-749.90681 

-2503.05607 

3 

-544.02252 

-1649.54460 

4 

-394.58328 

-1112.74704 

5 

-286.12668 

-775.13626 

6 

-207.42501 

-562.80681 

7 

-150.32551 

-429.26193 

8 

-108.90823 

-335.50205 

9 

-78.87451 

-262.01548 

10 

-57.10299 

-204.43719 

11 

-41.32715 

-159.34467 

12 

-29.90121 

-124.04846 

13 

-21.63011 

-96.43506 

14 

-15.64610 

-74.85272 

15 

-11.31906 

-57.99242 

16 

-8.19147 

-44.84499 

17 

-5.93108 

-34.60124 

18 

-4.29666 

-26.62805 

19 

-3.11304 

-20.44467 

20 

-2.25301 

-15.64648 

21 

-1.62421 

-11.94312 

22 

-1.15958 

-9.08705 

23 

-.81673 

-6.89740 

24 

-.56674 

-5.21708 

25 

-.38898 

-3.94327 

26 

-.26502 

-2.97364 

27 

-.17974 

-2.24154 

28 

-.12149 

-1.69079 

29 

-.08171 

-1.27761 

30 

-.05433 

-.96725 

31 

-.03517 

-.73306 

32 

-.02142 

-.55312 

33 

-.01123 

-.41080 

34 

-.00330 

-.29403 

35 

.00330 

-.19140 

36 

.00929 

-.11662 

37 

.01534 

-.06997 

38 

.02204 

-.04046 

39 

.03003 

-.01571 

40 

.03998 

0.00000 

41 

.05260 

.01571 

42 

.06871 

.04046 

43 

.08916 

.06997 

44 

.11478 

.11662 

45 

.14619 

.19140 

46 

.18359 

.29403 
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TABLE 1. - CONTINUED 





TABLE 1. - CONTINUED 


INDEX 

X 

z 

94 

163.52559 


95 

227.68509 


96 

317.04773 


97 

441.48577 


98 

614.73612 


99 

855.91313 



The LTRAN2-HI default grid consists of 113 x 97 points in the x,z plane 
and is listed in Table 2 and illustrated in Fig. 3. For the flat plate case 
the physical extent of the grid is +200c in x and +2327c in z with 48 grid 
points on the airfoil. The far-field boundary conditions in the program are 
also reflecting so the physical extent is quite large. The computational 5,n 
grid in LTRAN2-HI is related to physical x,z space by 


X = C 


z 


7 ^ 


The X grid is again fixed while the z grid depends only upon the thickness 
ratio. 


TABLE 2. LTRAN2-HI DEFAULT GRID FOR 6 = 0.005 


INDEX 

X 

Z 

1 

-200.00000 

-2326.50510 

2 

-132.08390 

-1694.05485 

3 

-87.24290 

-1240.86615 

4 

-57.65660 

-914.37169 

5 

-38.15274 

-677.87267 

6 

-25.31051 

-505.62532 

7 

-16.86743 

-379.48372 

8 

-11.32718 

-286.59564 

9 

-7.70010 

-217.81429 

10 

-5.33166 

-166.59850 

11 

-3.78896 

-128.24781 

12 

-2.78570 

-99.36814 

13 

-2.13252 

-77.49680 

14 

-1.70418 

-60.83830 

15 

-1.41785 

-48.07719 

16 

-1.21875 

-38.24505 

17 

-1.07500 

-30.62553 
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TABLE 2. - CONTINUED. 


INDEX 

X 

18 

-.95000 

19 

-.82500 

20 

-.70000 

21 

-.57500 

22 

-.45625 

23 

-.35000 

24 

-.25625 

25 

-.18125 

26 

-.12500 

27 

-.08188 

28 

-.05375 

29 

-.03625 

30 

-.02375 

31 

-.01500 

32 

-.00813 

33 

-.00250 

34 

.00250 

35 

.00750 

36 

.01250 

37 

.01750 

38 

.02250 

39 

.02750 

40 

.03250 

41 

.03813 

42 

.04563 

43 

.05500 

44 

.06500 

45 

.07500 

46 

.08563 

47 

.09875 

48 

.11703 

49 

.14204 

50 

.17188 

51 

.20313 

52 

.23438 

53 

.26563 

54 

.29688 

55 

.32813 

56 

.35938 

57 

.39063 

58 

.42188 

59 

.45313 

60 

.48438 

61 

.51563 

62 

.54688 

63 

.57813 

64 

.60938 

65 

.64063 


-24.68612 

-20.02908 

-16.35590 

-13.44147 

-11.11518 

-9.24715 

-7.73799 

-6.51133 

-5.50815 

-4.68267 

-3.99918 

-3.42973 

-2.95229 

-2.54946 

-2.20742 

-1.91512 

-1.66372 

-1.44609 

-1.25646 

-1.09014 

-.94330 

-.81278 

-.69600 

-.59080 

-.49539 

-.40827 

-.32815 

-.25398 

-.18482 

-.11988 

-.05848 

0.00000 

.05848 

.11988 

.18482 

.25398 

.32815 

.40827 

.49539 

.59080 

.69600 

.81278 

.94330 

1.09014 

1.25646 

1.44609 

1.66372 

1.91512 
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TABLE 2. - CONTINUED 
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The grid Implemented In XTRAN2L Is designed to Improve the accurecy of 
the calculated forces for two-dimensional flow problems. It Is 80 x 61 points 
in x,z and is listed in Table 3 and shown In Flq. 4. The physical arid 
extends ?20c In x and ?2Fc In z with 51 grid points lying on l^hr alrfoi 

^’’S'^XTRAN “'by ^ computltlojli 


X = c 

z = 




The computational K grid Is fixed for all problems while the n orld will vary 
depending upon M«, 6 and y*. The XTRAN2L grid takes advantage of the 
far-fl eld radiation boundary conditions incorporated in the program by 
reducing the physical extent covered. Reducing the area covered by the grid 
allows the number of points In the grid to be decreased, permits the inclusion 
vTolm^ points on the airfoil, and reduces the computational cost. The 
XTRAN2L grid consists of a smooth stretching In z and In x upstream and 
downstream of the airfoil. 

On the airfoil, the LTRAN2-NLR and LTRAN2-HI proarams concentrate grid 
points near the leading and trailing edges because' of rapid changes In 

Ulfrih® spacing on the airfoil is maximum 

^ nriddle of the airfoil. Yet the chanae 1n pressure across a shock is 
Of the same magnitude as the change In pressure near the airfoil 
edges so the grid spacing around the shock should also be fine. 
Unfortunately the shock very likely will occur In the region where both qrlds 
have relatively large grid spacing, causing reduced shock resolution. 

*u pressure gradient across a shock may be greater than that near 

the the XTRAN2L grid has no special concentration of points 

near the airfoil trailing edge. The XTRAN2L grid points on the airfoil were 
equally spaced to allow a better shock definition for ary position and 
transient motion. Keyfitz, Mel nick and Grossman^^ show that the grid 
spacing near a blunt leading edge Is important In TSP codes and the best 
results may be obtained for a medium mesh spacing. Thus no fine arid spacing 
Is used near the leading edge but a grid point was added near the' nose to aid 
in defining a blunt airfoil. 


TABLE 3. XTRAN2L DEFAULT GRID 


INDEX 

X 

z 

1 

-20.00000 

-25.00000 

2 

-16.30961 

-23.36111 

3 

-13.10054 

-21.77778 

4 

-10.33987 

-20.25000 

5 

-7.99453 

-18.77778 

6 

-6.03136 

-17.36111 

7 

-4.41705 

-16.00000 
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TABLE 3 . CONTINUED 


JDEX 

X 

Z 

8 

- 3.11817 

- 14.69444 

9 

- 2.10110 

- 13.44444 

10 

- 1.33204 

- 12.25000 

11 

-.77699 

- 11.11111 

12 

-.40170 

- 10.02778 

13 

-.17160 

- 9.00000 

14 

-.05175 

- 8.02778 

15 

-.00667 

- 7.11111 

16 

.00667 

- 6.25000 

17 

.02000 

- 5.44444 

18 

.04000 

- 4.69444 

19 

.06000 

- 4.00000 

20 

.08000 

- 3.36111 

21 

.10000 

- 2.77778 

22 

.12000 

- 2.25000 

23 

.14000 

- 1.77778 

24 

.16000 

- 1.36111 

25 

.18000 

- 1.00000 

26 

.20000 

-.69444 

27 

.22000 

-.44444 

28 

.24000 


29 

.26000 

-.11111 

30 

.28000 

-.02778 

31 

.30000 

0.00000 

32 

.32000 

.02778 

33 


.11111 

34 

.36000 

.25000 

35 

.38000 

.44444 

36 

.40000 

.69444 

37 

.42000 

1.00000 

38 

.44000 

1.36111 

39 

.46000 

1.77778 

40 

.48000 

2.25000 

41 

.50000 

2.77778 

42 

.52000 

3.36111 

43 


4.00000 

44 

.56000 

4.69444 

45 


5.44444 

46 

.60000 

6.25000 

47 

.62000 

7.11111 

48 


8.02778 

49 

.66000 

9.00000 

50 

.68000 

10.02778 

51 

.70000 

11.11111 

52 

.72000 

12.25000 

53 

.74000 

13.44444 

54 

.76000 

14.69444 

55 

.78000 

16.00000 


56 .80000 17.36111 













TABLE 3. CONTINUED 


INDEX 

>< 

z 

57 

.82000 

18.77778 

58 

.84000 

20.25000 

59 

.86000 

21.77778 

60 

.88000 

23.36111 

61 

.90000 

25.00000 

62 

.92000 


63 

.94000 


64 

.96000 


65 

.98000 


66 

1.00000 


67 

1.02000 


68 

1.12274 


69 

1.35473 


70 

1.75323 


71 

2.35080 


72 

3.17695 


73 

4.25898 


74 

5.62249 


75 

7.29170 


76 

9.28970 


77 

11.63860 


78 

14.35968 


79 

17.47352 


80 

21.00000 



The XTRAN3S program is a 3-dimensional TSP finite difference code. For 
analysis with the XTRAN2L program the default x,z grid defined at the root 
chord is studied. The grid coordinates are listed in Table 4 and the near- 
field grid is shown in Fig. 5. The physical extent of the grid is from -15c 
to 27c in X and +13c in z with 39 grid points on the airfoil. The XTRAN3S 
program has reflecting far-field boundary conditions and the grid is limited 
to 60 X 40 points in the x,z plane. This limitation results in the boundaries 
being defined very close to the airfoil as compared to the previously 
described grids. The physical grid is invariant and related to the XTRAN3S 
computational grid by 


x-Xle = ^ 
z = n 

where x.^ the leading edge location (=0 at root chord). The 

computatfffnal and physical grids at the root chord are identical and are 
problem independent. For analysis with XTRAN2L a grid row at z = 0 was added 
to the XTRAN3S default z grid because the XTRAN2L program requires a grid row 
at z = 0 while the XTRAN3S program does not. 
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TABLE 4. XTRAN3S DEFAULT GRID 


INDEX 

X 

Z 

1 

-15.37500 

-13.03750 

2 

-7.69250 

-6.63750 

3 

-3.85000 

-3.43750 

4 

-1.92750 

-1.83750 

5 

-.96500 

-1.03750 

6 

-.48250 

-.63750 

7 

-.24000 

-.43750 

8 

-.11750 

-.33750 

9 

-.05500 

-.28750 

10 

-.02250 

-.26250 

11 

-.00500 

-.23750 

12 

.00500 

-.21250 

13 

.01500 

-.18750 

14 

.02500 

-.16250 

15 

.03500 

-.13750 

16 

.04500 

-.11250 

17 

.06000 

-.08750 

18 

.08000 

-.06250 

19 

.10000 

-.03750 

20 

.13000 

-.01250 

21 

.16000 

0.00000 

22 

.19000 

.01250 

23 

.22000 

.03750 

24 

.25000 

.06250 

25 

.28000 

.08750 

26 

.31000 

.11250 

27 

.34000 

.13750 

28 

.37000 

.16250 

29 

.40000 

.18750 

30 

.43000 

.21250 

31 

.46000 

.23750 

32 

.49000 

.26250 

33 

.52000 

.28750 

34 

.55000 

.33750 

35 

.58000 

I .43750 

36 

.61000 

.63750 

37 

.64000 

1.03750 

38 

.67000 

1.83750 

39 

.70000 

3.43750 

40 

.73000 

6.63750 

41 

.76000 

13.03750 

42 

.79000 


43 

.82000 


44 

.85000 


45 

.88000 


46 

.91000 


47 

.93500 


48 

.96000 


49 

.98000 
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TABLE 4. - CONTINUED 


INDEX 

X 

z 

50 

1.00000 


51 

1.02500 


52 

1.07500 


53 

1.17500 


54 

1.37500 


55 

1.77500 


56 

2.57500 


57 

4.17500 


58 

7.37500 


59 

13.77500 


60 

26.57500 



A revised XTRAN3S x,z grid was designed to improve the accuracy of the 
calculated forces and is listed in Table 5 and shown in Fig. 6. The physical 
extent of the grid is +20c in x and +25c in z with 39 points lying on the 
airfoil. The grid also is problem-independent and has a z = 0 grid line added 
for the two-dimensional analysis. The grid was optimized for the available 
size of 60 x 40 points in the x,z plane. To minimize the reflections from the 
z boundary, the extent covered by the z grid is increased to +25 
chordlengths. A smooth stretching is used in z and in x upstream and 
downstream of the airfoil. The points on the airfoil are equi-spaced with an 
additional point near the leading edge for better nose definition. 

TABLE 5. REVISED XTRAN3S GRID 


INDEX 

X 

z 

1 

-20.00000 

-25.00000 

2 

-14.70824 

-22.50164 

3 

-10.47170 

-20.13478 

4 

-7.16272 

-17.89941 

5 

-4.65677 

-15.79553 

6 

-2.83281 

-13.82314 

7 

-1.57362 

-11.98225 

8 

-.76633 

-10.27285 

9 

-.30308 

-8.69494 

10 

-.08199 

-7.24852 

11 

-.00877 

-5.93360 

12 

.00877 

-4.75016 

13 

.02632 

-3.69822 

14 

.05263 

-2.77778 

15 

.07895 

-1.98882 

16 

.10526 

-1.33136 

17 

.13158 

-.80539 

18 

.15789 

-.41091 

19 

.18421 

-.14793 

20 

.21053 

-.01644 

21 

.23684 

0.00000 
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TABLE 5. - CONTINUED 


INDEX 

X 

z 

22 

.26316 

.01644 

23 

.28947 

.14793 

24 

.31579 

.41091 

25 

.34211 

.80539 

26 

.36842 

1.33136 

27 

.39474 

1.98882 

28 

.42105 

2.77778 

29 

.44737 

3.69822 

30 

.47368 

4.75016 

31 

.50000 

5.93360 

32 

.52632 

7.24852 

33 

.55263 

8.69494 

34 

.57895 

10.27285 

35 

.60526 

11.98225 

36 

.63158 

13.82314 

37 

.65789 

15.79553 

38 

.68421 

17.89941 

39 

.71053 

20.13478 

40 

.73684 

22.50164 

41 

.76316 

25.00000 

42 

.78947 


44 

.84211 


45 

.86842 


46 

.89474 


47 

.92105 


48 

.94737 


49 

.97368 


50 

1.00000 


51 

1.02632 


52 

1.19383 


53 

1.62332 


54 

2.42771 


55 

3.71531 


56 

5.59119 


57 

8.15791 


58 

11.51599 


59 

15.76424 


60 

21.00000 



RESULTS AND DISCUSSION 


Linear Flat Plate Analysis 
LTRAN2-NLR GRID 

The first grid analyzed using the pulse-transfer function technique Is 
the LTRAN2-NLR grid. The grid was analyzed usinq the reflecting far-fleld 
boundary conditions used In the LTRAN2-NLR program. Figure 7 prkents ct 
versus reduced frequency for this grid at M» = 0.850 and for the 
corresponding exact values obtained from a kernel function solution. In 
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addition to obtaininq the unsteady forces for the complete TSP equation on 
this (jrid, the unsteady forces are also calculated for the low frequency 
equation and boundary conditions'^ and are included to show the effect of 
adding the time derivative terms to the equations. 

Roth finite difference solutions exhibit spurious oscillations in the 
calculated unsteady forces. The complete equation solution exhibits slightly 
larger oscillations than the low frequency results. The low frequency 
solution deviates from the complete and exact solutions for reduced 
frequencies greater than 0.075. Several modifications of the grid were 
conducted to determine the cause of this oscillation. The resulting unsteady 
forces are very insensitive to the x grid extent and spacing off of the 

airfoil while the z grid extent and spacing are very critical. For the 
physical time span covered in the flat plate cases, the maximum distance from 
which disturbances can be reflected back to the airfoil is 41c in x and 77c in 
z. The x grid boundaries were moved to within 4 chordlengths of the airfoil 
before any noticeable change occured in the unsteady forces, while changes 

were noticeable when the z boundaries came within 75 chordlengths of the 
airfoil . 

The oscillations are due to internal numerical reflections in the grid 
caused by stretching too rapidly in the z direction and were traced to the 
grid spacing at approximately 9 chordlengths from the airfoil. The 

disturbances from the airfoil due to the pulse are unable to propagate 

accurately between the grid points which are approximately 2 chordlengths 

apart and are partially reflected back to the airfoil. The reflected 
disturbances cause oscillations in the force time histories as shown in Fig. 1 
and result in the oscillations in the force frequency response functions. 

Adding the time derivative terms by going from the low frequency to the 

complete equation causes the oscillations to increase in magnitude. 

As the reduced frequency increases, the LTRAN2-NLR forces begin to 
deviate from the exact values. This is not due to the grid but is caused by 
using too large a time step. The calculated forces can be brought into closer 
agreement with the kernel function values at the higher reduced frequencies by 
decreasing the time step. To maintain accuracy at the low reduced 

frequencies, the number of steps should be increased such that the total time 
the solution spans remains constant. The phenomenon of decreasing agreement 
with the exact kernel function values for increasing reduced frequency for the 
pulse-generated forces occurs for all cases in this paper since they all were 
run for the identical number of time steps and total time. 

LTRAN2-HI 


Results from the LTRAN2-HI grid analyzed using reflecting far-field 
boundary conditions are shown in Fig. 8. This grid is very similar to the 
LTRAN2-NLR grid and shows the same oscillations in unsteady forces. Internal 
grid reflections affect the forces for reduced frequencies above 0.20. 
The magnitude and extent of the oscillations in the forces are not as severe 
for the LTRAN2-HI grid as for the LTRAN2-NLR grid. This is due to the 

slightly finer spacing of the LTRAN2-HI grid around 9 chordlengths away in z 
than that found in the LTRAN2-NLR grid. 


17 



XTRAN2L Grid 


The results for the XTRAN2L grid are shown in Figs. 9 and 10. In order 
to see the effect of the radiation boundary conditions, the grid was first 
analyzed with the reflecting far-field boundary conditions used for the 
LTRAN2-NLR and LTRAN2-HI grids (Fig. 9). The calculated forces agree quite 
well with the kernel function values except for the oscillations near 
k = 0.08. These oscillations are due to reflections of disturbance from the z 
boundary at +25 chordlengths from the wing. Including the radiation far-field 
boundary conditions alleviates the oscillations as shown in Fig. 10. 

In addition to giving more accurate forces, the XTRAN2L grid gives a 
considerable savings in computational cost due to the fewer number of grid 
points used. The computer time required to calculate forces using this grid 
is 34% less than the time taken for the LTRAN2-NLR grid and 53% less than that 
for the LTRAN2-HI grid. The program size is also reduced by 15% and 28% as 
compared to LTRAN2-NLR and LTRAN2-HI programs respectively. 

XTRAN3S Default Grid 


The forces calculated for the default XTRAN3S root chord x,z grid using 
reflecting far-field boundary conditions are shown in Fig. 11. Very poor 
agreement with the exact values is exhibited. The XTRAN3S grid exhibits both 
internal and boundary reflections. For k < 0.3 the oscillations are due to 
reflections from the z boundary and for k > 0.3 are caused by internal grid 
reflections from approximately 3 chordlengths above and below the airfoil. 

XTRAN3S Revised Grid 


The results for the revised XTRAN3S root chord x,z grid using reflecting 
far-field boundary conditions are shown in Fig. 12. The revised grid exhibits 
much better agreement with the exact values than the original grid. The 
amplitude and frequency of the boundary reflections are decreased due to the 
larger extent of the grid in z. The effects of internal grid reflections are 
eliminated for k < 0.50 and minimized for k > 0.50. Neither the boundary nor 
the internal reflections can be eliminated completely due to the program size 
restrictions and reflecting boundary conditions. The results represent a 
compromise grid designed to optimize accuracy for k < 0.50, which is the range 
of interest for flutter analysis. 

6% Parabolic Arc Analysis 


To investigate non-linear transonic effects, calculations were made for a 
six percent thick parabolic arc airfoil using the LTRAN2-NLR, LTRAN2-HI, and 
XTRAN2L grids. The LTRAN2-NLR and LTRAN2-HI grids were analyzed using 
reflecting far-field boundary conditions while the XTRAN2L grid used the 
radiation boundary conditions. The forces calculated using the XTRAN2L grid 
at M«, = 0.850 are shown in Fig. 13.: To verify the accuracy of the 

pulse-transfer function technique for this non-linear case, harmonic 
oscillatory solutions were obtained for discrete reduced frequencies. Figure 
13 demonstrates that the harmonic and pulse unsteady forces agree well with 
each other for moderate reduced frequencies and verifies that for this case, 
the unsteady forces may be treated as locally linear. As demonstrated in the 
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linear flat plate cases, the agreement becomes worse with increasing reduced 
frequency. The agreement can he improved by decreasino the time step for the 
pulse case. 

The forces for the three grids are compared in Figs. 14 and 15. Figure 
14 demonstrates that all three grids predict values of C£ that are close 
to one another. The LTRAN2-NLR and LTRAN2-HI grids cause internal reflections 
as in the linear flat plate case. The magnitude of the oscillations for this 
non-linear case are smaller than for the corresponding linear cases shown in 
Figs. 7 and 8. The pitchino moments shown in Fig'. 15 also exhibit the 
internal reflections and, in addition, demonstrate a difference in the value 
of C|n^ between the three grids. 

To illustrate the source of the difference exhibited in the values of 
Cma» the steady-state pressure distributions of the airfoil for each grid 
are shown in Fig. 16. As the grid spacing on the airfoil becomes finer, the 
shock definition becomes better and the shock strength increases. As the 
shock strength increases, the lift is not changed to any measureable degree 
but the moment increases. Grid spacing along the airfoil can thus be very 
important to the shock definition and hence the resulting forces. 


CONCLUDING REMARKS 

A pulse-transfer function method for calculatinq unsteady aerodynamic 
forces for a wide range of reduced frequencies has been implemented in an 
unsteady transonic small perturbation finite difference code. The forces were 
determined for a two-dimensional linear flat plate and compared to exact 
theoretical values to evaluate the grid. The LTRAN2-NLR grid exhibited inter- 
nal reflections which can be eliminated by improving the z stretching. The 
LTRAN2-HI grid was shown to be superior to the LTRAN2-NLR arid but to still 
have internal reflections. The addition of far-field radiation boundary 
conditions for the complete equation permitted a decrease in arid size and 
resulted in a 3 ^% savings in computing time while eliminating the internal 
reflections. The unsteady forces calculated with the XTRAN3S arid were shown 
to be inaccurate due to boundary and internal numerical reflections. A 
redesigned grid utilizing the same total number of points but spanning a 
larger extent and having an improved stretching in z was shown to increase the 
accuracy. 

Calculations for a 6% thick parabolic arc airfoil demonstrated that arid 
problems encountered for the linear flat plate airfoil also occur for 
non-linear cases. The pulse-transfer function technique was shown to work for 
a non-linear case, indicating that the unsteady forces can be treated as 
locally linear. In addition, the grid spacing was shown to influence the 

shock definition which in turn affects the unsteady pitching moment. Finer 
grid spacing on the airfoil was used to adequately capture shock definition 
and motion. 
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Fig. 2 LTRAN2-NLR default grid near the airfoil for M« * 0.850, 6 = 0.005. 
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Fig. 3 LTRAN2-HI default grid near the airfoil for 6 = 0.005. 
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Fig. 4 XTRAN2L default grid near the airfoil. 
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Fig. 5 XTRAN3S default root chord x,z grid near the airfoil 
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Fig. 6 Revised XTRAN3S root chord x,z grid near the airfoil. 
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Fig. 7 UnsteacI/ forces calculated for a flat plate airfoil using the 
LTRAN-NLR default grid, M« = 0.850. 
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Fig. 9 Unsteady forces calculated for a flat plate airfoil using the XTRAN2L 
default grid with reflecting far-fleld boundary conditions. 
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REDUCED FREQUENCY, k 

Fig. 10 Unsteady forces calculated for a flat plate airfoil using the 

XTRAN2L default grid with radiation far-field boundary conditions, 
Me = 0.850. 
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Fig. 11 Unsteacly forces calculated for a flat plate airfoil using the 
default JCTRAN3S root chord x,z grid, M« = 0.850. 
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REDUCED FREQUENCY, k 

Fiq. 12 Unsteady forces calculated for a flat plate airfoil usinq the 
revised XTRAN3S root chord x,z grid, = 0.850. 
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Fig. 15 CjUgj for a 6 % parabolic arc airfoil calculated using the 
LTRAN2-NLR, LTRAN2-HI and XTRAN2L default grids, M« = 0 
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Fig. 16 Steady-state pressure distributions for a 6% parabolic arc airfoil 
calculated using the LTRAN2-NLR, LTRAN2-HI and XTRAN2L default 
grids, M» = 0.850. 
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